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Abstract 



Quasi-saddles or inherent saddles of the potential energy surface, U, of a liq- 
uid are defined as configurations which correspond to absolute minima of the 
pseudo-potential surface, W = |Vt/| 2 , as identified by a multi-dimensional min- 
imisation procedure. The sensitivity of statistical properties of inherent saddles 
to the convergence criteria of the minimisation procedure is investigated using, 
as a test system, a simple liquid bound by a quadratically shifted Lennard- Jones 
pair potential with continuous zeroth, first and second derivatives at the cut- 
off distance. The variation in statistical properties of saddles is studied over a 
range of error tolerances spanning five orders of magnitude. The largest value 
of the tolerance corresponds to that used for the unshifted LJ liquids in a pre- 
vious work (J. Chcm. Phys. 115, 8784 (2001)). Based on our results, it is 
clear that there are no qualitative changes in statistical properties of saddles 
over this range of error tolerances and even the quantitative changes are small. 
The lowest magnitude eigenvalue, |wg|, of the Hessian is, however, found to be 
very sensitive to the tolerance; as the tolerance is decreased, |u>q| is found to 
show an overall decrease. This indicates that if convergence criteria are not 
strict, absolute or low-lying minima of W(r) will be diagnosed as having no 
inflexion directions. The results also show that it is not possible to set up an 
unambiguous numerical criterion to further classify the quasi-saddles into true 
saddles which contain no zero curvature, non-translational normal modes and 
inflexion points which have one or more zero-curvature normal mode directions. 
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1 Introduction 

A number of theoretical approaches have been developed to understand the 
relationship between the potential energy surface (PES) and the properties of 
condensed phases. Since the potential energy surface, U(r), is a complicated 
function of the 3N dimensional position coordinates of an TV-atom system, most 
approaches have concentrated on identifying features of U (r) most relevant to 
understanding a particular phenomenon. The most widely used of such energy 
landscape approaches is inherent structure analysis which focuses on the sam- 
pling of local minima of the PES H§, ![§. To perform the inherent structure 
sampling, configurations are sampled from a well-defined ensemble and are re- 
ferred to as instantaneous configurations in this work. A steepest-descent or 
equivalent local minimisation procedure is then initiated for each instantaneous 
configuration in order to obtain the nearest local minimum. All instantaneous 
configurations which quench to the same minimum or inherent structure are 
referred to as belonging to the basin of that minimum. In general, the phase 
with lowest entropy and internal energy will correspond to location of the sys- 
tem in the basin of the global minima and phase transitions will correspond 
to occupancy of higher energy metastable minima. Inherent structure analysis 
has proved to be an extremely useful tool in understanding melting, the glass 
transition and protein folding |4|, ||, [SJ, |7[ 0, [| [uj[ O]. It is, however, limited 
by the fact that it focuses entirely on minima whereas an understanding of dy- 
namical phenomena requires a knowledge of factors which govern the passage 
of the system between basins of different minima. For example, the dynamical 
slowing down associated with the glass transition requires an understanding of 
why rate constants for inter-basin transitions become essentially zero on ex- 
perimental time scales. In this context, it has recently been suggested that 
saddle configurations of the PES, which mark the border between adjoining 
basins of minima, must play a crucial role. Such saddle configurations will have 
Vf = and will show negative curvature in one or more orthogonal direc- 
tions. Saddles of order 1 with just a single negative curvature direction will 
correspond to a simple transition state connecting two adjacent minima. These 
stationary points of U(r) will also correspond to absolute minima of the pseudo- 
potential surface, W(r) = i|VC/| 2 fT|, g 0, fj]| [u], ■ To identify the saddle 
points which play a significant role in the system dynamics, it has been sug- 
gested that inherent saddles, analogous to inherent structures, can be identified 
by a steepest descent minimisation from an instantaneous configuration to the 
nearest local minimum on the pseudo-potential surface, W(r)= i|VC/| 2 , where 
U(r) is the potential energy as a function of the 3N-dimensional position vec- 
tor r. The statistics of inherent saddles has been recently analysed for several 
model glass- formers, including a binary Lennard- Jones mixture(BLJ), a modi- 
fied monoatomic Lennard- Jones(MLJ) liquid and the binary soft-sphere (BSS) 
mixture. Simulations show that these systems undergo a transition from saddle- 
dominated to minima-dominated dynamical regimes at the mode-coupling tran- 
sition temperature, Tmc- For T > TmCi the system is localised largely in border 
regions and basin-hopping is facile. Below Tmc, the average saddle order goes 
to zero implying that the system occupies interiors of the basins of local minima 
and basin-hopping becomes an activated process. The saddle order is found to 
be an essentially linear function of the saddle configurational energy. By ex- 
trapolating the linear trend one can obtain a threshold energy, Uthr, at which 
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the saddle order goes to zero. The average configurational energy at the mode- 
coupling temperature is found to be very close to Vthr- Since a glass may be 
viewed as a very high- viscosity liquid, it is expected that many of the statistical 
features of the stationary points of the PES of glass-formers will be shared by 
simple liquids. Recent studies of the stationary points of Lennard- Jones and 
Morse liquids verify this conjecture jl8| [ljj. The statistical properties of inher- 
ent saddles are shown to carry interesting signatures stcming from variations in 
generic features of the interatomic potential, such as the range and curvature; 
for example, the configurational energy is shown to depend linearly on saddle 
order with a slope that is proportional to the range of the pair potential. These 
recent studies therefore provide strong evidence that stationary points of the 
potential energy surface constitute a physically significant set of topographic 
features of the energy landscape. 

The computational strategy for identifying stationary points or inherent sad- 
dles in the above studies has been to locate the minima of the pseudo-potential 
surface, W(r) = ^|VJ7| 2 ; within the limits of computational precision, the min- 
ima for which W(r) = are taken to be stationary points p(| . As pointed out 
in ref . , this approach must, however, be applied very carefully because any 
minimisation algorithm will sample not only absolute minima of W(r) which 
correspond to the stationary points of the true potential U(r) but also low- lying 
minima of W(r) which correspond to inflexion points of U(r) At an in- 

flexion point, the Hessian matrix, H, must have one or more zero eigenvalues; 
the eigenvectors corresponding to these zero eigenvalues do not correspond to 
basin-crossing displacements. The question therefore arises as to whether the 
physical significance of inherent saddle analysis is compromised because, given 
finite numerical precision limits, it may be difficult to distinguish between true 
saddles and inflexion points. The available numerical studies of this issue are 
limited but they indicate that the number of inflexion directions is always very 
small and does not significantly affect the statistical estimates of saddle prop- 
erties [|l6| [I?]]. For example, in the case of the 256 atom MLJ system studied 
in ref.14, no point with more that four inflexion directions was found. In the 
liquid regime, the ratio of inflexion directions to negative curvature directions 
was approximately 0.1. It rises to 0.8 very close to Tmc but even then no major 
shift in the estimated value of Tmc was seen. 

Given the usefulness of stationary point analysis in understanding liquid- 
state dynamics, we feel it is worthwhile to perform a detailed computational 
study of a model system. In this work, we therefore investigate the sensitivity 
of statistical properties of inherent saddles to details of the minimisation pro- 
cedure. For this purpose, we define a simple liquid bound by a quadratically 
shifted Lennard- Jones pair potential with continuous zeroth, first and second 
derivatives at the cut-off distance. This ensures that no loss of precision in the 
minimisation procedure is present due to discontinuities in the pair potential 
at a finite cut-off distance. Previous studies of stationary points of bulk liquids 
and glasses have used cither unshifted pair potentials with long-range correc- 
tions |l^, [l9| or quadratically shifted pair potentials with continuous zeroth 
and first derivatives but discontinuous second derivatives jl^, [l3|, [l4|, [l5|, |l~6| ]. 
The only exception in the literature known to us is the simulation of the bi- 
nary soft-sphere system but this may represent a special case because of the 
absence of any attractive tail in the pair potential p7| . Unlike in other stud- 
ies, we consider a one-component bulk liquid, rather than a binary mixture; 
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these two classes of systems differ only in the ease with which crystal nuclcation 
takes place but the crucial statistical features of stationary points are common 
to both glass-formers and simple liquids. This monoatomic liquid bound by a 
quadratically shifted Lennard- Jones potential (QSLJ) defines a test system for 
which we can set reasonably demanding accuracy limits for the minimisation 
procedure. The low-lying minima of the pseudo-potential, W(r) = i|Vl/| 2 , as 
located by the minimisation procedure are referred to as quasi-saddles; however, 
unless the context requires distinguishing clearly between true saddles and in- 
flexion points, we also refer to quasi-saddlcs as saddle configurations. Changes 
in various properties of the saddle configurations are analysed as a function of 
the precision limits set for the minimisation algorithm. The possibility of setting 
up an unambiguous numerical criterion to distinguish between true saddles and 
inflexion points is also discussed This allows us to obtain a useful benchmark 
of the degree of error introduced by a given level of convergence in the min- 
imisation procedure. The results should be useful in two contexts. Firstly, it 
allows one to estimate the computational effort to be expended in minimisation 
given the saddle properties of interest. Secondly, in a number of contexts, it is 
of interest to study liquids which do not have a finite range e.g. the Lennard- 
Jones and other well-known model systems. If the statistical analysis of saddle 
configurations is to be a generally applicable tool to study the disordered phase 
(liquids and glasses), it is essential to establish to what extent such an analy- 
sis will be limited by computational problems associated with the minimisation 
procedure. The paper is organised as follows. Section 2 contains details regard- 
ing the computational methods employed. Results are presented in Section 3. 
Section 4 contains the conclusions. 



2 Computational Methods 
2.1 Potential Energy Surface 

The potential energy surface (PES), U(r), is defined as a scalar function of the 
3N-dimensional position vector, r. The PES studied in this work is pair-additive 
i.e. 

N 

»=1 3<i 

where is the distance between atoms i and j and the pair potential, it p (r), 
is a quadratically shifted Lennard-Jones function of the form 

u p (r) = 4e |(-) _ (-) j+ar 2 + 6r + c for r < R c 

= for r > R c . (2) 

The a, b and c parameters are determined by requiring that u p , du p /dr and 
d 2 Up/dr 2 be zero at r = R c . We take e, a and m as reduced units for energy, 
length and mass respectively; all our results are given in reduced units. The 
pair potential defined in equation (2) is referred to as the QSLJ potential in this 
paper. The two-parameter quadratically shifted Lennard-Jones potential used 
in previous studies of saddles of binary LJ mixtures corresponds to a shifting 
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function ar 2 + c and is referred to as the SLJ potential. The monoatomic SLJ 
system was originally discussed in ref. (22| . We fix the spherical cut-off distance 
at 2.5cr when computing the shifting parameters. Figure 1 compares the QSLJ, 
SLJ and unshifted Lennard-Jones (LJ) potentials. 



2.2 NPT Ensemble Monte Carlo Simulations 

All the simulations reported in this paper were performed using the isothermal- 
isobaric ensemble Monte Carlo method of McDonald |2j| . The partition function 
of an iV-particle system in this ensemble may be written as: 



n pp 

WNPT 



A 3N Nl 



J d(lnv)v N+1 exp(-/3Pv) J dsexp(-(3V(s;L)) (3) 



where P is the external pressure, /3 = l/fc^T is the inverse temperature, A = 
yjh 2 /2nmkBT is the thermal de Broglie wavelength, v and L are the volume 
and length of the simulation cell respectively and s = r/L is the 3N-component 
position vector scaled with respect to the length of the simulation cell. The 
volume, v, and the particle positions, r, are the Metropolis variables. The 
random walk is performed using two types of moves: (i) volume moves which 
sample in hif and (ii) particle moves which displace the spatial coordinates 
of a single, randomly chosen particle. In the case of the unshifted Lennard- 
Jones system, long-range corrections were employed and the spherical cut-off 
distance, R Cl was scaled whenever a volume move was made; the average value 
of R c ss 2.8cr. For the shifted potentials, the cut-off distance was fixed at 2.5<7. 

A simulation cell with 125 particles and rhombic dodecahedral boundary 
conditions was used in the Monte Carlo simulations. The solid phase structure 
taken to be a face-centred cubic (fee) lattice. Volume moves constituted 7.4% 
of trial moves. Acceptance ratios for volume and particle moves were kept at 
approximately 50% . Run lengths were between two and four million configura- 
tions with cquilibriation periods of one million to two million. This was found to 
be sufficient to converge configurational averages, such as the average potential 
energy, density and bond orientational parameters, to better than 1%. In order 
to obtain statistical averages over inherent saddle configurations, 100 instan- 
taneous configurations were sampled from the NPT-MC simulations at equal 
intervals and used to generate inherent saddles. The simulations were carried 
out at a reduced pressure of P* = 0.67. At this pressure, the melting point of 
the unshifted Lennard-Jones system has been established from previous work 
as T* = 0.75 §1. 



2.3 Locating Inherent Saddles 

Locating inherent saddles requires one to find the absolute minima for the 
pseudo-potential function, W(r) = ^\VU\ 2 . Any local minimisation technique 
equivalent to a steepest descent procedure can be used to perform the min- 
imisation. In a previous study | fl8| , three minimisation procedures for locating 
inherent saddles were compared: (i) the conjugate gradient algorithm given in 
ref. [^5| . (ii) the variable metric LBFGS algorithm, as implemented by Liu and 
Nocedal [|26| and (iii) Powell's direction set method which requires no gradient 
information. The LBFGS algorithm was found to be computationally the most 
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efficient and was therefore used in this study to locate local minima in W(r) 
using VW = — H • F. The termination condition for the LBFGS algorithm is 

<rmax(l, \\ Xi \\) (4) 

where gi and Xj are the gradient and solution vectors for the i-th iteration and 
1 1 • 1 1 denotes the Euclidean norm. Clearly the more demanding the accuracy 
requirements i.e. the smaller the error tolerance t, the more computationally 
demanding the minimisation procedure will be. This is illustrated in Table 1 
which shows the average number of iterations required to converge to the nearest 
saddle for a given value of r for three different temperatures. The error tolerance 
r is the key convergence parameter with respect to which we wish to test the 
sensitivity of properties of quasi-saddles. 

While absolute minima of W(r) must correspond to true stationary points 
of U(r), local minima in W(r) need not correspond to true stationary points of 
Z7(r) with all elements of the gradient vector, VU, equal to zero. Some local 
minima may be inflexion points where both the first and second derivatives in 
some directions are zero. In principle, such flex points can be distinguished 
from true saddles by looking at the absolute magnitude of W(r). However, 
given finite numerical precision it may be difficult to unequivocally distinguish 
between low-lying and absolute minima. Alternatively, such inflexion points 
can be located by examining whether any of the eigenvalues of the Hessian 
for a configuration corresponding to a minimum in W(r), other than the three 
translational modes, is unusually small. However, the lowest magnitude, non- 
translational, eigenvalue of the Hessian is very sensitive to the accuracy of the 
minimisation procedure. Relaxed tolerances can lead to dramatic undercounting 
of the number of flex points while rigorous convergence criteria lead to almost 
99% of local minima of W(r) being identified as inflexion points (for example, 
one can contrast the results presented in refs.]l2| and fig] . 



3 Results and Discussion 

3.1 Isobars of the Quadratically Shifted Lennard- Jones 
liquid 

The behaviour of the qudratically shifted Lennard- Jones system (QSLJ), as de- 
fined in Section 2.1, is mapped out along the P* = 0.67 isobar to facilitate 
comparison with the unshifted Lennard-Jones system. We have also performed 
NPT ensemble simulations at the same pressure for the two-parameter quadrat- 
ically shifted Lennard- Jones (SLJ) potential with discontinuous second deriva- 
tive, as defined by Stoddard and Fordp2|. The binary Lennard-Jones mixtures 
for which inherent saddle analysis has been performed recently employ the SLJ 
type of cut-off conditions. Figures 2(a) and 2(b) compare the temperature de- 
pendence of the configurational energy per particle and the number density 
respectively along the P* = 0.67 isobar for the QSLJ, SLJ and LJ systems. 
The thermodynamic melting temperature, T* n , is known to be 0.75 for the LJ 
system at this pressure [[24J . The metastability limit of the solid from our NPT 
simulations is = 0.88; such a 15-20% difference between T* % and T£ is fairly 
typical. Though the difference in the well depth of the pair potential between 
the LJ and SLJ systems is very small (sec Figure 1), the T£ values for the SLJ 
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system is about 13% lower than that of the unshifted LJ system because of the 
finite-range of the SLJ pair potential. The QSLJ system has an even lower value 
of Tj* = 0.64 because of the reduced well-depth. At any given temperature, the 
QSLJ system has the lowest number density and highest configurational po- 
tential energy. Since we are interested in the stationary points sampled in the 
disordered liquid phase, the results of the inherent saddle analysis discussed 
below refer to liquid state simulations for T* > 0.75. 

3.2 Sensitivity of Saddle Properties to Accuracy of Min- 
imisation Procedure 

A saddle point is characterised by its configurational energy, U s , and its index 
density, n s , which corresponds to the fraction of imaginary modes; the NPT 
averages are denoted by (U) s and (n) s . In the work of Broderix et al, (U) s and 
(n) s are referred to as parametric averages with respect to T as the parameter. 
For the binary LJ mixture, they showed that the points ((U) s , (n) s ) lie on a 
straight line. They also defined the geometric average, U s (n s ) as the average 
configurational energy of all saddles of index density n s and showed that the 
straight line fits to the geometric and parametric averages are identical, within 
the limits of statistical error. This linear relationship between the saddle order 
and configurational energy has been subsequently shown to be true for Lennard- 
Jones and Morse liquids; moreover, the slope, dU s /dn s has been found to be 
proportional to the range of the potential |l^, [lj|. The linear form of the 
(U) s ((n) s ) curves in the disordered, liquid phase can be understood on the 
basis of a localised mechanism for generating imaginary frequency displacement 
modes which requires rearrangement of only a few neighbouring atoms. 

Figure 3(a) shows the parameteric points ({U) s , (n) s ) obtained from minimi- 
sations corresponding to three different values of the error tolerance, r = 10~ 10 , 
10 -12 and 10~ 15 . It is well know that error tolerances for systems with finite 
cut-offs and discontinuous scroth, first or seond derivatives cannot be made very 
small. The tolerance value of e = 10 -10 corresponds to that used in a previ- 
ous study of the unshifted LJ liquid |jg|| . The tolerances used for the Morse 
potentials in a more recent study are similar. The slopes of the straight line 
fits are almost unchanged over five orders of magnitude in r. The results at 
r = 10~ 12 and 10~ 15 are indistinguishable. For t = 10~ 10 , for a given value 
of the saddle energy, the saddle order is somewhat overestimated. Figure 3(b) 
shows the temperature dependence of the ensemble averaged index density. The 
behaviour with t is similar to that seen for the U s (n s ) curves. Figure 3(c) shows 
the ensemble-averaged saddle configurational energies as a function of temper- 
ature. The convergence with r is much better for the saddle energies than for 
index densities. The reason is that the convergence of the near-zero eigenvalues 
of the Hessian is relatively strongly affected by the error tolerance compared to 
the overall saddle configurational energies. To substantiate this point, we have 
compared the change in the absolute magnitude of the non-translational eigen- 
values, \ujf\, of the Hessian with r for several configurations. Figure 4 shows the 
values of the lowest, \u>q\, and third lowest, magnitude eigenvalues of the 
Hessian for an arbitrary sequence of configurations for r = 10~ 14 and 10~ 15 . 
The lowest magnitude eigenvalue varies by over a factor of 10 when t is reduced 
by an order of magnitude. There is, on average, a decrease in the magnitude of 
lUq as t decreases. The magnitudes of the third-lowest eigenvalue (and all other 
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higher magnitude eigenvalues) are practically unchanged with r. The sensitiv- 
ity of the lowest magnitude eiegenvalues of the Hessian to the error tolerance 
of minimisation is responsible for the relatively large error in (n) s compared to 
(U) s discussed above. 

3.3 Distinguishing between True Saddles and Inflexion Points 

We have approached the possibility of setting up an unambiguous numerical 
criterion to separate the quasi-saddles into true saddles and inflexion points by 
comparing the location of each quasi-saddle on a two-dimensional plot with x— 
and y— axes corresponding to W(r s ) and \loq\Jljo\ respectively where luq is the 
lowest magnitude non-translational eigenvalue of the Hessian (see Figure 5) and 
is the Einstein frequency. True saddles should correspond to very small 
values of W and relatively large values of |u>o| and therefore should cluster in 
the upper left hand corner of the graph. Inflexion points should correspond to 
relatively larger values of W and near-zero values of |wq| and therefore should 
cluster in the lower right hand corner of the graph. If two such well-defined 
clusters of points are seen, then an unambiguous numerical distinction between 
true saddles and inflexion points can be made. Figure 5 shows, however, that 
the situation is quite the opposite. Results for two different temperatures are 
shown and are representative of the behaviour in other simulation runs. Only 
a few points (less than 5 in 100) can be classed unambiguously as true saddllcs; 
the others form an almost continuous distribution. This is presumably why in 
ref.16 the use of tight convergence criteria for minimisation led to an identifica- 
tion of 99% of the quasi-saddles as inflexion points. We suggest that the reason 
for this behaviour is that in the absence of symmetry constraints requiring some 
collective normal mode directions of quasi-saddles to have zero curvature, the 
distinction between saddles and inflexion points is not physically very mean- 
ingful because it attempts to distinguish between zero and near-zero curvature 
situations. Coupled with the fact that the number of near-zero curvature modes 
is always small, as can be seen even from the saddle normal mode spectra de- 
scribed in ref. [jl8] as well as from the results in ref. [|l6| , it would appear that the 
effect of reducing t will be to introduce small corrections in the saddle index 
density. Since the physical significance of such quasi-saddles stems from the 
fact that the majority of negative curvature directions will correspond to basin- 
crossing displacements, improving the error tolerance will not qualitatively alter 
any results. 

4 Conclusions 

The purpose of this study was to test the sensitivity of statistical properties 
of quasi-saddles of the potential energy surface, U(r), to the convergence cri- 
teria used for minimisation on the pseudo-potential surface, W(r) = \WU\ 2 . 
For this purpose a quadratically shifted Lennard- Jones (QSLJ) pair potential 
was defined which has continuous zcroth, first and second derivatives at the 
cut-off distance. Inherent saddles (or quasi-saddles) of the bulk QSLJ liquid 
were studied as a function of the error tolerance, r, of the LBFGS minimisation 
algorithm. The variation in statistical properties of saddles was studied over 
a range of r values covering five orders of magnitude. The largest r values 
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correspond to those used for unshifted LJ and Morse liquids simulations previ- 
ously. Based on our results, it is clear that there are no qualitative changes in 
statistical proeprties of saddles over this range of error tolerance and even the 
quantitative changes are small. The ensemble-averaged saddle index densities 
were found to be more sensitive to r than the corresponding averages over saddle 
configurational energies. The lowest magnitude eigenvalue, |wo| 2 , of the Hessian 
was, however, found to be very sensitive to the tolerance; as r was decreased, 
|wo| 2 , was found to show an overall decrease. This indicates that if convergence 
criteria are not strict, absolute or low- lying minima of W^r) willl be diagnosed 
as having no inflexion directions, as was found in ref.12 and 18. For the smallest 
value of r, we also considered the possibility of setting up a numerical criterion 
to clearly distinguish between true saddles and inflexion points by considering 
the value of both W(r) as well as |w 2 | at a quasi-saddle identified by minimisa- 
tion on the pseudo-potential surface. We find however, that the quasi-saddles 
do not fall into two such clear and distinct categories. The reason appears to 
be that of the 3N normal mode directions associated with a low-lying or abso- 
lute minima of W(r), a few directions will correspond to near-zero curvatures. 
Deducing which of these near-zero curvatures is actually identically equal to 
zero and therefore corresponds to an inflexion direction is necessary in order 
to distinguish between true saddles and inflexion points Moreover, this appears 
not to be a physically very meaningful distinction since the physically impor- 
tant feature of quasi-saddles is the existence of negative curvature directions 
indicating their location at the border between adjacent minima. 
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Table 1 

Average Number of iterations of the LBFGS minimisation algorithm for a 
given error tolerance, r, at a reduced temperature T* for the quadratically 
shifted Lennard-Jones (QSLJ) liquid. 
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54 


1077 


10895 


1.29 


35 


898 


17362 
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Figure Captions 



1. Comparison of the quadratically shifted Lennard- Jones (QSLJ) pair po- 
tential with continuous zeroth, first and second derivatives with the two- 
parameter quadratically shifted Lcnnard-Jones (SLJ) potential defined in 
rcf.22 as well with the unshifted Lennard-Jones (LJ) potential. The well- 
depth, e, and size parameter, a, of the LJ potential are taken as reduced 
units of energy and length respectively. 

2. Dependence of configurational properties on reduced temperature T* at 
a reduced pressure P*=0.67 for the bulk LJ, SLJ and QSLJ systems: (a) 
average potential energy per particle < U* > and (b) reduced density p* 

3. Statistical properties of saddle configurations for three different values of 
the error tolerance, t, used in the LBFGS minimisation algorithm: (a) 
Parametric averages of saddle configurational energies, (U) s , versus index 
densities, (n) s (b) saddle index densities as a function of the reduced tem- 
perature, T* and (c) saddle configurational energies, (U) s , as a function 
of the reduced temperature, T*. Configurational energies U are in units 
of e per particle. 

4. Sensitivity of the lowest magnitude, \u>q\, and third lowest magnitude, 
eigenvalues of the Hessian for a set of 20 saddle configurations sampled 
from the simulation at T* = 1.05 for for r = 10~ 14 and 10~ 15 . The y-axis 
has a logarithmic scale. 

5. Scatter plots showing the correlation of the |wq|/cj|. values with the value 
of the pseudo-potential, W(r) = |VC/| 2 , for saddle configurations sampled 
from simulations at two different temperatures where loe is the Einstein 
frequency. 
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